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Abstract 

We perform the Hamiltonian analysis for the lowest-order effective action, 
up to second order in derivatives, of the complete Hofava theory. The model 
includes the invariant terms that depend on di IniV proposed by Bias, Pujolas and 
Sibiryakov. We show that the algebra of constraints closes. The "Hamiltonian" 
constraint is of second-class behavior and it can be regarded as an elliptic partial 
differential equation for N. The linearized version of this equation is a Poisson 
equation for N that can be solved consistently. The preservation in time of the 
Hamiltonian constraint yields an equation that can be consistently solved for a 
Lagrange multiplier of the theory. The model has six propagating degrees of 
freedom in the phase space, corresponding to three even physical modes. When 
compared with the \R model studied by us in a previous paper, it lacks two 
second-class constraints, which leads to the extra even mode. 
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1 Introduction 



There has been a debate about the consistency of the Hofava theory pQ, which is a 
proposal for a UV completion of general relativity that could be a renormalizable theory. 
Some of the papers involved in this debate can be found from Refs. [2] to [22]. The 
discussion has been mainly focused on the presence of a physical mode additional to 
the ones of general relativity that could be the potential source of instabilities. This 
extra mode is present at all scales in the projectable version of the theory since in this 
case there is no local Hamiltonian constraint (see [2] for the Hamiltonian analysis on 
the corresponding second-order actionjE On the other hand, we found jl] that the 
lowest-order truncation of the original Hofava theory is physically equivalent to general 
relativity for arbitrary values of the constant A. This result encourages us to deepen 
the study on the nonprojectable formulation of the Hofava theory, to which this paper 
is devoted. 

Among the papers focused on the nonprojectable case, in Ref. [5] perturbative com- 
putations signaled the problem of the strong coupling of the theory (in particular of 
the extra mode), which implies the breakdown of its low-energy perturbative expansion. 
This behavior was confirmed in Ref. [6J. However, the extra mode did not manifest 
itself (as a propagating mode) in the linear-order perturbative analyses of Refs. [7, 8^. 
This apparent contradiction is explained by the argument given in [5J, where it was 
indicated that the extra mode is excited in perturbative analyses (at linear order) only 
on time-dependent and spatially nonhomogeneous background^]. 

In Refs. [TTJ [121 [131 CHI SI IS] the constraints of the original nonprojectable theory 
have been studied, as well as its lower-derivative truncations. The authors of Ref. [TT] 
met several difficulties when performing an early Hamiltonian analysis. They could not 
close the algebra of constraints and suggested the arising of an excessive number of con- 
straints (a result that has some similarity with [H]). Actually, these inconveniences are 
circumvented once the Hamiltonian constraint is regarded as a second-class constraint 
with no associated gauge symmetry [T3J HJ [15]. Nevertheless, the authors of Ref. [11] 
noticed the fact that the phase space of the original Hofava is of odd dimensionality. 
In [12] it was proposed that the nonprojectable theory is in general inconsistent due 
to the presence of a constraint that is absent in general relativity. However, for the 
models studied in Refs. [U [15] this constraint can be solved in closed ways. Still on 

In the projectable case the extra mode could be eliminated by enhancing the gauge symmetry 
group, see [5] 

4 In Ref. [7] some "deformed" versions of the theory were analyzed; but such deformations, although 
representing violations of the detailed balance principle, are just redefinitions of some constants in the 
potential. We prefer to keep their denomination as the "original" Hofava theory, in the sense that their 
potentials are written in terms of only the spatial curvature tensors (and the cosmological constant). 
In general we consider the several coupling constants as independent parameters; hence we do not use 
the detailed balance principle. 

5 In addition, in Ref. [9] the physical behavior of the extra mode in the projectable formulation of 
the theory has been analyzed and other modifications/extensions have also been considered (see, for 
example, [TO]). 
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Hamiltonian grounds, in Ref . [13] , assuming standard asymptotic behavior for the grav- 
itational variables (the same assumptions we are going to use in this paper), it was 
argued that the lapse function must be vanishing at least asymptotically in the general 
theory and everywhere for the lowest-order truncation. The same posture was assumed 
in Ref. |14j. Regarding this, in Refs. [U [TS] we performed Hamiltonian analysis for 
low-energy effective models for the Hofava theory, starting in [I] with the lowest-order 
model, 



As we mentioned above, in that paper we found that the model (11 .11) is physically 
equivalent to general relativity regardless of the value of the constant A, including the 
fact that it propagates two physical modes. This is due to the fact that, when performing 
Dirac's algorithm for the preservation of constraints, the condition K = emerges as 
a second-class constraint of the model (instead of having a vanishing lapse function, as 
was proposed in [13]). Thus, the term XK 2 drops out and the theory resulting from 
(II. ip evaluated on the constrained phase space coincides with general relativity in the 
particular gauge K = 0. Next, in Ref. [15] we incorporated an R 2 term in the potential. 
We found that the algebra of constraints closes, but some peculiarities arise. First, 
the R model lacks one of the second-class constraints of the model which leads 

to a physical odd-dimensional (five dimensions) submanifold in phase space, as was 
previously indicated in Refs. [TTJ [6]. This means that there is a physical mode with 
a first-order evolution equation, hence propagating half of the usual Cauchy data. It 
seems that this half mode will be persistent when adding higher-order (z = 3) terms in 
the potential. This is so because the lacking of one second-class constraint of (11.11) is 
a consequence of the explicit dependence on the lapse function of the constraint that 
would generate it. Whenever a (second-class) constraint depends on the lapse function, 
it is to be expected that its preservation in time generates an equation for a Lagrange 
multiplier of the theory, rather than a new constraint. This feature does not arise in 
( II. ip . but probably will hold in general in the original Hofava formulation with higher- 
order terms. Second, N is determined by a first-order partial differential equation (PDE) 
of flow type. The boundary data compatible with such an operator must be given on 
a noncharacteristic surface defined by a certain vector field. In any case, the equation 
is compatible with the expected asymptotic behavior of the variables that yields the 
Minkowski space at infinity. Moreover, in the perturbative analysis around a weakly 
varying background we did in Ref. [15], we found that the extra mode decouples at the 
lowest orders, hence recovering the physics of the model ( II. ip smoothly. The results of 
Refs. [U [15] were corroborated by the perturbative Hamiltonian analysis of Ref. |16j . 
As we mentioned above, the presence of the extra mode in Hofava theory was ana- 

6 When one considers N and its conjugated momentum as part of the canonical variables, the con- 
straints of the model Ijl.ip are the momentum constraint and four second-class constraints: the vanishing 
of the momentum of N, the Hamiltonian constraint, the vanishing of K, and an elliptic equation for 
N. The last three are the (ordered) chain of constraints yielded by imposing the preservation in time 
of previous constraints. These constraints reduce from 14 to four the number of canonical variables. 




(1.1) 
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lyzed in Ref . [6] . Those authors first made a counting of degrees of freedom in a model 
having up to a fourth-order term (R^R^) in the potential. The counting was formulated 
as a problem of initial Cauchy data in the corresponding equations of motion. Second, 
they focused on the physical behavior of the extra mode, but truncating to a second- 
order model with the aim of simplifying the computations. From their analysis they 
concluded that the extra mode suffers from very fast instabilities and strong coupling 
at the deep IR. However, it turns out that the second-order model they used to analyze 
the propagation of the extra mode is exactly the model (11. ip . which does not contain 
the extra mode (this was already pointed out by us in Ref. [15J). Based upon this and 
the smooth decoupling found in [T5J , we suggest being cautious with the behavior of the 
extra mode of the nonprojectable Hofava theory. 

The authors of [6] suggested that the inconveniences they found with the extra 
mode could be cured by the inclusion of additional terms depending on the spatial vector 
dj = di In N. These terms were not considered by Hofava in his original formulation. The 
proposal was formally presented by them in Ref. [TTj. Several terms that are invariant 
under the action of foliation-preserving diffeomorphisms, which are the underlying gauge 
symmetries of the Hofava theory, can be formed with appropriate combinations of a«, 
its derivatives and curvature tensors. In particular, the expressions Vja* and a^a 1 are 
second-order terms that are invariant under the foliation-preserving diffeomorphisms. 

The physical consistency of the Hofava theory with the terms of Bias, Pujolas and 
Sibiryakov has also been the subject of study. The problem of the strong coupling was 
reported to persist in Ref. [18]. However, a response from Bias, Pujolas and Sibiryakov 
[T9] (see also [20] ) pointed out that the complete theory can be made free from the strong 
coupling problem if the energy scale of its UV physics (where higher- derivative terms 
become important) is below the Planck scale. In Ref. [2T] it was confirmed that the 
theory exhibits strong coupling if no new scales are introduced. Perturbative analyses 
of this theory with cosmological interest have been made in Ref. [22]. 

Independently of the presence or not of undesirable physics behavior for the extra 
mode, specially taking into account that it decouples smoothly at the largest distances 
in the original Hofava model [15], one may be in favor of including the aj terms by 
following an argument of completeness of the theory. When dealing with effective the- 
ories it is normally assumed that all the terms of the order in consideration that are 
compatible with the gauge symmetries must be included in the Lagrangian. Therefore, 
if the quantum consistency of the Hofava could be proved, then the terms should be 
included since they will arise in quantum corrections^. 

With the aim of gaining more insight into the canonical structure of the nonpro- 
jectable Hofava theory, in this paper we study the Hamiltonian formulation of the 
lowest-order effective action (second-order in derivatives) that includes the Via* and 
did 1 terms. A previous Hamiltonian analysis of the nonprojectable Hofava theory with 
the terms of Bias, Pujolas and Sibiryakov was presented in Ref. [23]. In this paper we 
further develop the analysis of the constraints. In Sec. II we shall find all the constraints 

7 The renormalizability of the Hofava theory has been studied in Ref. [23] 



4 



of the theory explicitly and we shall see that Dirac's algorithm for extracting the con- 
straints of the theory ends with an elliptic PDE for a Lagrange multiplier of the theory. 
This closure follows the same behavior we have mentioned: the last constraint being 
generated is the one that depends on the lapse function. This analysis will allows us 
to count the number of physical degrees of freedom of the model. In Sec. Ill we shall 
show an interesting canonical transformation that is defined in terms of a conformal 
transformation and that leads to some simplifications in the canonical action. 

2 Analysis of the constraints 

The theory is written in terms of the Arnowitt-Deser-Misner variables 

ds\ = (-N 2 + NiN^dt 2 + 2N i dx i dt + g ij dx i dx j . (2.1) 

We denote by 7r u the momentum conjugated to gij and by <fi the one of N. After 
performing the Legendre transformation it turns out that the vector Ni can be regarded 
as a Lagrange multiplier, which is to be expected because of the gauge symmetries of 
the theory. This is not the case for N, which, together with its conjugated momentum, 
must be considered as part of the canonical variables. Hence the unconstrained phase 
space is parametrized by {(gij,^), {N, </>)}. tt^ and (f) behave as scalar densities under 
spatial coordinate transformations. 

We want to study the model as a local theory of gravity without topological effects. 
Therefore we assume that the whole spatial manifold is noncompact. In our varia- 
tional calculus all the configurations are fixed at spatial infinity and, since we do not 
include a cosmological constant, correspond to Minkowski space-time. We also assume 
that the field variables have the same asymptotic behavior of general relativity [25]: in 
asymptotically flat coordinates, as r approaches to infinity, they behave as 

9ij = 5 tl + 0(r- 1 ), 7r« = 0(r- 2 ), 

(2.2) 

N = l + 0(r^), Ni = 0{r~ l ). 

Note that for the consistency of the canonical action it is assumed that and N go as 
0{r~ 2 ). Indeed, this is the asymptotic behavior of g^ that is consistent with the one of 
7T 1 - 7 dictated by (12.21) . In particular, these conditions typically lead to finite functionals 
of the field variables, as, for example, the kinetic term J d 3 x7r lJ g^. 
The object 

ai = dilnN (2.3) 

transforms as a vector under the spatial sector of the foliation-preserving diffeomor- 
phisms and as a scalar under the time transformations. That is, under Diff(.M, F) 
transformations defined by St = f(t) and 5x l = ( l (t,x), 0,4 transforms according to 

5a { = C 3 djdi + di^cij + fhi . (2.4) 
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Therefore, it can be consistently used jTT] to add terms to the original Lagrangian of 
the Hof ava theory. 

The second-order action (without a cosmological constant), which is the most general 
one preserving time and spatial parity, is given by 



S = J dttPx. y /gN(G ijkl K ij K kl + R + aa^a' + /3V;a* 



(2.5) 



where a and are arbitrary coupling constants and 



K* = ^(^•-2V ( ^. ) ), (2.6) 

G ijM = \(g ik g jl + g il g jk )-\g ij g hl . (2.7) 
For the case of A ^ 1/3 the inverse of G^ kl is given by 

Gijkl = -j[9ik9ji + gugjk) — 3^ _ ^ 9ij9ki ■ (2-8) 

Asymptotically (I2.2p . y/gNVid 1 is of order 0(r~ 3 ) whereas ^/gNaia 1 is 0(r~ A ). 
Therefore, although one can express locally the former as the negative of the latter 
plus a spatial divergence, their contributions to the potential are definitely different due 
to the finite contribution of y/gNVia 1 at the spatial infinity. After the integration by 
parts, we get the final form of the Lagrangian, 



J d z x^gN{G iikl K ij K kl + R + aa^) + f3§ N 



(2.9) 



$ w = & dZidiN , a = a-(3. (2.10) 



S = I (If 

where 



To get this form of $at we have used the asymptotics (12.21) . Since the flux $jy is a surface 
term at spatial infinity, its functional derivatives vanish provided the variations have an 
asymptotic behavior of 0(r~ 2 ), the one consistent with the asymptotic assumptions in 
(12. 2p and on and N. Consequently, there are no contributions of $at to the field 
equations or to the Poisson brackets. In fact, Poisson brackets are distributions acting 
on functions of compact support; hence the Dirac delta and any of its derivatives acting 
on such functions at the boundary give a vanishing contribution. 

We obtain the Hamiltonian from the action (12.91) . The momentum 7r u has the 
universal form 

?l = G ijkl K kl . (2.11) 

This implies that the velocities gij can be completely solved in favor of 7r y if the con- 
dition A 7^ 1/3 is verified. We assume this condition throughout this paper. After 



6 



performing the Legendre transformation and adding the primary constraint 0, we get 
the Hamiltonian 

H = J d 3 x(Nn + NiW + acf)) - (3(P N , 

TT^TT kl _ (2 12) 

Ti = Qijki — — y/g{R + aaia 1 ) , v ' ' 

V 9 

TV = -2Vj7T j \ 

where Qijki is the inverse of G^ kl and a is a Lagrange multiplier that transforms as a 
scalar under spatial coordinate transformations. 

The Legendre transformation automatically incorporates the primary constraint TV 
into the Hamiltonian. However, notice that TV generates spatial coordinate transfor- 
mations only in and . In order to carry out a clear treatment, we would like 
to compute brackets with the full generator, which must include also the generator of 
the spatial coordinate transformations on iV and <fi (in any case, we can always go to 
the gauge Ni = and forget about these generators, but we prefer to perform a gen- 
eral treatment). The generator of infinitesimal spatial coordinate transformations on N 
and (j) is (pdtN, which vanishes in the constrained surface. Therefore, we redefine the 
momentum constraint by 

U* = H i + (j)d l N (2.13) 

and replace H l by 1-L 1 in the Hamiltonian 02. 12j) . A similar consideration was included 
in Ref. [21]. 

To ensure the preservation in time of the <fi constraint we need to compute only its 
Poisson bracket with J cPxN'H. We obtain 

{(f), H} = -n- 2a-yfg(V i a i + a^) . (2. 14) 

Therefore, the density 

n = n + 2a^(V i a i + aid) (2.15) 

is a secondary constraint of the theory. Again, the expression ^fg~N{Via l + a,ia l ) is 
an exact divergence whose integral gives the flux of at infinity. Thus, by solving 
"H in terms of % and replacing it into the Hamiltonian, which adds another boundary 
term proportional to $tv, we may write the Hamiltonian as a sum of constraints and a 
boundary term. After doing this we arrive at the final form of the Hamiltonian 

H = J d 3 x(NH + NiW + ct0) - 0<f> N , (2.16) 

H = g^i^^ + Vgi-R + ^Via' + aia')), (2.17) 
yd 

= 2a -(3, (2.18) 
and ?V is given in (12.131) . 
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If we use variations in the strong sense as in [23] , Hamiltonian (" 12 . 16[) is functionally 
differentiable if and only if = 2a, that is, j3 — 0. The remaining boundary term is 
then relevant in this mathematical sense. From the physical point of view, it contributes 
to the gravitational mass of the theory. Thus we see that, under the variational scheme 
of Ref. [23], the V^a* term must be excluded from the Lagrangian. However, we notice 
that under the asymptotic conditions we are considering, where cjij and N behave as 
0(r~ 2 ), the arbitrary variations Sgij must be of the same order (D(r~ 2 ). Otherwise the 
integrability of the Lagrangian is violated. Assuming then this behavior, it turns out 
that the variation of the action gives rise to a boundary integral in terms of the variations 
which become zero. The functional differentiability of the action is then satisfied without 
adding boundary terms a la Regge and Teitelboim. In any case, we think it would be 
the quantum theory which will determine if there are quantum contributions or not 
to the boundary term. The coefficient beta will then be determined from quantum 
corrections. An important point related to it is the quantum stability of the action, a 
problem directly related to the positiveness of the gravitational mass. This problem has 
been recently addressed in Ref. J2JJ]. The quantum stability of the theory will determine 
bounds to the coefficient beta. 

Condition % = is a nonlinear elliptic PDE for TV; hence in principle we can solve 
it for N. Indeed, for the variable a/A this condition becomes a linear, homogeneous, 
elliptic PDE, 

(ij kl\ 
45V 2 -R + g ijk i^-J VN = . (2.19) 

To check the compatibility with the boundary conditions, at the asymptotic limit we 
expand vA = 1 + n and cast Eq. (12.191) as an equation for n. We get that the dominant 
terms are 

4adidin = R. (2.20) 

Since R is of 0(r~ 3 ), this equation dictates n = (9(r _1 ), which is the assumed behavior 
in O. 

It is worth studying a bit more the existence and uniqueness properties of Eq. (12.191) . 
since A is a dynamical variable of the theory which we expect to be fixed by the Hamil- 
tonian constraint. Any source of indetermination on N could lead to inconsistencies 
of the theory, or may force us to reinterpret the Hamiltonian constraint as a condition 
for another variable. The existence and uniqueness of the solution of Eq. (I2.19p are 
rather nontrivial aspects because of the presence of the R and tt 2 terms in the operator. 
We may give a partial answer to this question by linearizing the equation for all field 
variables, which can be achieved by a perturbative analysis around a Minkowski back- 
ground. We expand \J~N — 1 + n, gij = 8ij + hij and ir l i is considered of first order in 
perturbation. Under these settings, a perturbative expansion around Minkowski space- 
time is equivalent to going to the asymptotic limit. Hence we may read the linearized 
version of Eq. (I2.19P easily from Eq. (I2.20p . obtaining 

AadidiH = didjhij — didih . (2.21) 
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This is an ordinary Poisson equation in flat space, subject to the Dirichlet boundary 
condition = and with asymptotic decay of the source of 0(r~ 3 ), which is fast 
enough to ensure the finiteness of the solution. Therefore, the solution of Eq. (I2.2ip 
with the given boundary condition exists and is unique. 

Moreover, since the differential operator is elliptic, we can give a result on nonper- 
turbative grounds about the existence and uniqueness of the solution of the PDE (I2.19p . 
If the nonderivative terms of Eq. (I2.19p satisfy the condition 

sgn(a) (R - g- l g ijkl ^n kl ) > (2.22) 

in all the spatial submanifold, then the Lax-Milgram theorem can be used to obtain 
directly that the weak solution of ( 12 . 1 9[) subject to the prescribed boundary condition 
exists and is uniqu^l It would be nice to prove the existence and uniqueness of the 
weak solution when the condition (I2.22p is relaxed. 

From the previous experience [Ij with the model without the cij terms, one expects 
that and "H are of second-class behavior. Hence it is useful to have the Poisson brackets 
between them. The nonzero brackets, evaluated on the constrained phase space, are 
equal to 

{ J d 3 xeH, J d 3 y V H} = j d 3 x ^ ~ ^ (rjVh)* + 2a{vd i e)g ijkl a j ir kl ) - (e rj) , 

(2.23) 

{J d 3 xeH, J d 3 yri<i)} = 2a J d 3 x^gr]{V 2 t - V i (eo i ))JV _1 . (2.24) 

Now we demand the time preservation of H. This requires computing its Poisson 
brackets with j d 3 xN'H and J d 3 xa<j), which can easily be read from (I2.23P and f)2.24p . 
We obtain 



{U, H} = 2a^/g[V 2 (a/N) + a%(cr/N)} - 2aNg ijkl [V\aH kl ) + 2aV 



TT kl ] 



+^^iv- 1 v i (iv 2 W). (2 ' 25) 

Therefore, the preservation in time of % implies the following equation for a 

a[V 2 {a/N) + a'diia/N)) = J , (2.26) 

J = - (J^jj ^V 4 (iV 2 VV) + a^g tjkl [V\a^ kl ) + 2a l a^ kl ] . (2.27) 

Equation (12. 26ft is a second-order, linear, elliptic PDE for a with a source term 
independent of a. Thus we can solve this equation for a in a unique way if the operator 



8 With weak we mean a solution in the sense of distributions: if L stands for the operator acting on 
y/N, then there exists a unique function that satisfies the boundary conditions and J d 3 xfLy/~N = 
for any smooth function / of compact support . 
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is compatible with the prescribed boundary conditions. First, by varying the canonical 
action with respect to (f) we get the relation N = a. From this relation we may read the 
asymptotic behavior of a, but being careful with the subtleties arising with the time 
dependence. We recall that in (I2.2p it is assumed that asymptotically is of 0(r~ 2 ). 
By going back to Eq. ( I2.20p . we note that n is of 0(r~ 2 ). Therefore, we conclude that 
a = 0(r~ 2 ). Now we analyze the asymptotic behavior of Eq. (12.261) . In the source J the 
asymptotically leading term is the first term, which decays as 0(r~ A ). The dominant 
term for a in the left-hand side is oc didi<j\ hence we get 

adAa = £>(r~ 4 ) , (2.28) 

which is satisfied by a = 0(r~ 2 ). Since we can solve Eq. (I2.26P for the Lagrange 
multiplier a, Dirac's algorithm for the constraints ends at this step. 

We have ended up with the momentum constraint H 1 and the constraints 4> and H. 
From (I2.23P and (I2.24p we confirm that and % are the second-class constraints of the 
theory. The physical degrees of freedom are given by 

(# Physical D.O.F.) = 

(2.29) 

Canonical var.) - 2x (| 1 st class const.) — 2 nd class const.) = 6 . 

This leaves us with six independent degrees of freedom in the canonical space, which are 
equivalent to three propagating even modes. Two of them correspond to the graviton 
and the remaining one is an even scalar mode. Such a scalar mode is activated in the 
metric gr^, as was focused in Ref. [TT] . 

By going to the a = limit we get that Eqs. (I2.19P and (I2.26P becomes, respectively, 

-R + Qijki = 0, 2.30 

9 

(A — l)Vi(iV 2 VV) = 0. (2.31) 

We see that iV decouples from (I2.30P and a from (12.31 p . The first of these equations 
defines the Hamiltonian constraint of the XR model [3]. The only solution of the second 
equation that satisfies (I2.2p is n = 0, which is interpreted as a second-class constraint in 
this limit. The preservation in time of ir yields the second-class constraint (V 2 — R)N = 
0. Finally, the preservation of this last constraint yields an equation for a of the type 
(V 2 — R)a = J . This limit is physically equivalent to general relativity [1]. 



3 A conformal transformation depending on N 

One might be wondering about to what extent the spatial derivatives of iV arising in the 
action, Hamiltonian and constraints can be simplified by a conformal transformation of 
the metric depending on N. To focus this issue, we start by noting that the conformal 
transformation 

& = JV 2e <to, 7r« = JV-V, (3.1) 
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where e is an arbitrary real parameter, is a canonical transformation if we compensate 
it with the transformation 

= - 2eiV~ 1 7r . (3.2) 

Indeed, we get 

7r% + 0A> = #% + 0A>, (3.3) 
and the nonvanishing Poisson brackets are 

{~g tj (t,x),n kl (t,y)} = + 8$)8{2 - y) , 

(3-4) 

{iV(t,x),0(t,y)} = 5(x-y). 

These brackets were computed in terms of the old variables, but from now on they can 
be redefined in terms of the canonical variables (#y,7r u ) and (N,<j>), (13. 4p being the 
fundamental brackets. Now the primary constraint = takes the form 

^ = -2eiy _1 7f . (3.5) 

After performing the transformations (13. ip - (13. 2p . we get that the term J d 3 xN% 
of the Hamiltonian given in (I2.12p becomes 

J d 3 xNH = J d 3 x^N (^N 3e g ijk i^— - N~ e (R + x^^a^ - 4e$ w , (3.6) 

where \ = a + 2e 2 — 4e. Thus, we can drop the term with derivatives of N out of the 
bulk integral in (13.61) by putting \ = 0, which gives 



e = 1 ± VI - a/2 (3.7) 

under the condition 

a<2. (3.8) 

Therefore, for the sector of the space of parameters where (I3.8P is verified, the (bulk) 
Hamiltonian can be written without any explicit dependence on the derivatives of N. 
However, derivatives of iV still arise when one evaluates the condition for the preservation 
in time of the constraint (13. 5p . that is, the Hamiltonian constraint. Curiously, there is a 
special value of the coupling constant a for which the Hamiltonian gets an even simpler 
form. For a = 2 the formula (13. 7p yields e = 1 and (13 .6p takes the form 

J d 3 xNU = J d 3 xy^ (^N 4 g ijkl ^— -Rj- 4$ N . (3.9) 

We may see that in this case there is not dependence on N of the potential. This implies 
in particular that the Hamiltonian constraint acquires a totally algebraic dependence 
on N. This constraint becomes 

g m m kl N* = ~gR. (3.10) 

It would be interesting to elucidate the role of special values of the coupling constants, 
such as a = 2, from the point of view of the renormalization-group flow of the theory. 
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4 Conclusions 



We have performed the Hamiltonian analysis to the lowest-order effective action of the 
complete nonprojectable Hofava theory. Following the proposal of Bias, Pujolas and 
Sibiryakov, this action contains a term proportional to (<9jlniV) 2 and a boundary term 
which comes from the V 2 iV term. This model has the momentum constraint as the 
first-class constraint. The second-class constraints are the vanishing of the momentum 
conjugated to iV and the constraint needed for its preservation in time, which is anal- 
ogous to the so-called Hamiltonian constraint in general relativity. These constraints 
reduce to six the number of independent canonical variables, corresponding to three 
physical modes. It is remarkable that the extra scalar mode in this model is even, that 
is, evolves with a second-order equation. This fact was one of the motivations of Bias, 
Pujolas and Sibiryakov for introducing the extra terms in Ref. [T7]. This even scalar 
must be contrasted with the odd extra scalar found in Refs. 0, [TB] for models with 
quadratic-curvature terms in the potential and without Oj-dependent terms. 

When compared with the lowest-order effective action of the original Hofava theory, 
which is the XR model studied by us in Ref. [3|, the model studied in this paper lacks two 
second-class constraints. From the Hamiltonian-analysis point of view, this is the origin 
for the extra even scalar mode absent in the \R model, which is physically equivalent 
to general relativity. 

The Hamiltonian constraint is a second-order elliptic PDE for iV totally compatible 
with the standard (flat) asymptotic behavior of all the gravitational variables. In order 
to give support for the strategy of solving for N this constraint, we have linearized it 
and have found that it becomes a Poisson equation that can be solved in a closed way for 
N. Moreover, if the nonderivative terms satisfy a kind of non-negativity condition, the 
standard theory of PDEs ensures that the solution for N exists and is unique, at least in 
the sense of distributions. This is a totally nonperturbative result. It will be interesting 
to give the general proof of the existence and uniqueness of the solution for N. The 
preservation in time of the Hamiltonian constraint yields a second-order elliptic PDE 
for the Lagrange multiplier a of the theory compatible with the asymptotic conditions, 
hence ending Dirac's algorithm at this step. In this equation the operator acting on 
a does not contain nonderivative terms, which makes it simpler to conclude that the 
solution exists and is unique. 

We have also shown that for a particular value of one of the coupling constants the 
Hamiltonian constraint acquires an algebraic dependence on N after a suitable canonical 
transformation. 

Upon these results we conclude that the lowest-order truncation of the nonpro- 
jectable Hofava theory with the terms of Bias, Pujolas and Sibiryakov has a consistent 
and closed algebra of constraints, as was indicated in Ref. [24J, but the theory man- 
ifestly deviates from Einstein's general relativity at the deep IR, unlike the original 
Hofava theory [4J. For the viability of the theory beyond the mathematical consistency, 
the physics of the extra scalar mode should be contrasted with the phenomenology, as 
was concluded by Bias, Pujolas and Sibiryakov in Ref. |17j . 
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The consistent structure we have found for the Hamiltonian formulation of the effec- 
tive theory depends crucially on the form of the operators acting on N in the Hamiltonian 
constraint and on o in its corresponding equation. We have seen that these operators 
are of second-order elliptic form. However, from our analysis it cannot be asserted that 
the Hamiltonian formulation will remain consistent when higher-order operators are in- 
cluded, particularly the z = 3 terms required for renormalizability, since they raise the 
order of the mentioned operators. For example, this will happen in the Hamiltonian 
constraint when the fourth-order term (did 1 ) 2 is included. A promising feature about 
this point is that, if spatial parity is imposed as preserved symmetry, odd-dimensional 
operators are forbidden [T7]. Thus, the elliptic form of the equations for N and a 
possibly will be preserved by the higher-order operators. 

The model we have studied in this paper can be extended in several ways. An 
immediate step would be the inclusion of a cosmological constant in the Lagrangian. In 
this case the Hamiltonian analysis should be started from the very beginning since the 
boundary conditions and the asymptotic behavior of the field variables must be modified. 
We do not expect any modification in the number of constraints that emerge from Dirac's 
procedure, since the Hamiltonian constraint will depend on N and its preservation in 
time will generate an equation for a. Moreover, this equation will remain intact since a 
cosmological-constant term has no derivatives, hence it does not contribute to the source 
J of a. However, the details about the compatibility between these equations and the 
boundary conditions should be investigated carefully in order to arrive at a definitive 
conclusion about the consistency of the Hamiltonian formulation of the theory with 
cosmological constant. 

* * * 



Note added: While we were preparing the first version of our manuscript, the paper 
"Hamiltonian Structure of Hofava Gravity", by W. Donnelly and T. Jacobson [27], was 
published on the Web. With regard to the second-order action, that paper contains 
the Hamiltonian and the constraints, and we agree with their results. They also argued 
that the preservation in time of the Hamiltonian constraint leads to a condition for the 
Lagrange multiplier a (called Nw by them). We have found explicitly this condition in 
Eq. (I2.26p . In addition, they found that for the asymptotically flat case the Hamiltonian 
can be written as a sum of constraints plus a boundary contribution. We also found 
this structure for the Hamiltonian, but also the contribution of the V^a* term at the 
boundary. However, we must stress that, if variations in the strong sense of Ref. [25] 
are used, then the V^a* term cannot be included in the Lagrangian in order to obtain a 
functionally differentiable action with those variations. 
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